Coherence and pairing in a doped Mott insulator: Application to the cuprates 
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The issues of single particle coherence and its interplay with singlet pairing are studied within 
the slave boson gauge theory of a doped Mott insulator. Prior work by one of us^ showed that 
the coherence scale below which Landau quasiparticles emerge is parametrically lower than that 
identified in the slave boson mean field theory. Here we study the resulting new non-fermi liquid 
intermediate temperature regime characterized by a single particle scattering rate that is linear in 
temperature (T). In the presence of a d-wave pair amplitude this leads to a pseudogap state with 
T dependent Fermi arcs near the nodal direction. Implications for understanding the cuprates are 
discussed. 

PACS numbers; 



In a recent paper— we developed a unified phenomeno- 
logical picture of some of the most striking properties 
of the underdoped cuprates. A key aspect of the phe- 
nomenology is the existence of a low 'coherence' scale 
Tcoh for the single particle excitations below which Lan- 
dau quasiparticles become well defined. We took T^oh 
to rise as the doping x is increased (roughly order the 
superconducting transition temperature in zero mag- 
netic field), and to be only weakly affected by moderate 
field strengths. The purpose of this paper is to explore 
a specific microscopic model where the issue of electron 
coherence and its interplay with the singlet pairing (that 
undoubtedly also exists in the underdoped cuprates) can 
be studied in detail. 

It has long been recognized that the properties of the 
cuprates must be understood within the framework of 
doping a Mott insulator'^. A useful theoretical framework 
to discuss the doped Mott insulator is provided by a slave 
particle treatment of the "i — J" model. In the standard 
"slave boson" theory of the t — J modelS.;^, the mean field 
phase diagram looks as shown in Fig. [T^. The slave boson 
field h condenses at a temperature Ti, ^ Tc ^ x which 
is then identified with Tcoh- In addition d-wave singlet 
pairing of the spins sets in at a T* line that decreases 
with increasing doping. Fluctuations effects beyond the 
mean field theory have also been extensively addressed. 
At temperatures T > T^, non-fermi liquid physics has 
been found^i^. 

In this paper we revisit the slave particle gauge theory 
and show, following recent work by one of usi, that true 
single particle coherence is not established till an energy 
scale Tcoh that is parametrically lower than the mean 
field estimate T;,. The existence of an intermediate tem- 
perature non-fermi liquid regime in between and Tcoh 
was pointed out relatively recently^, and was missed in 
the prior literature. Ref. 'l! briefly described some prop- 
erties of this regime which we dub an "Incoherent Fermi 
Liquid" (IFL). Here we study it in greater detail in the 
specific context of the cuprates. In the absence of pair- 
ing (i.e above the T* line), the IFL is characterized by 
a T^/^ specific heat, a constant spin susceptibility, and a 
sharp Fermi surface leading to enhanced "2Kf" response 
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FIG. 1: Schematic slave boson phase diagram of the doped 
t — J model, (a) The mean field phase diagram (b) Phase 
diagram beyond mean field. All dashed lines represent 
crossovers. 



in both spin and charge channels. Most interestingly we 
show that the single particle spectrum possesses peaks 
at the Fermi surface with a scattering rate oc T unlike in 
a Landau fermi liquid. The crossover to the true Fermi 
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liquid happens only at Tcoh ^ x's 't^ Ti,. Below T* in the 
regime dubbed PG (for pseudogap) in Fig. [1)3, this scat- 
tering persists despite the development of a d-wave pair 
amplitude and leads to a pseudogap and T-dependent 
Fermi arc behavior. 

All the prior theoretical work'^ on the strange metal 
and pseudogap phases within the slave boson frame- 
work has focused on the non-fcrmi liquid (NFL) and spin 
gap(SG) regimes of Fig. [T]d. IFL or its descendant PG 
have not been studied thus far but may be more relevant 
to the temperatures actually probed in experiments. In- 
deed it was recognized early on^ that the slave boson 
mean field theory predicts a Tb (and hence a coherence 
temperature) that is far too high (order lOOOif or higher 
at optimal doping). Thus the lowering of the true co- 
herence scale described above from the slave boson mean 
field theory is potentially significant, and the IFL, PG 
regimes need to be explored. 

To describe our results, first consider the phase dia- 
gram in the absence of pairing, i.e above T*. Then the 
system can be viewed as resulting from doping a spin liq- 
uid Mott insulator with a Fermi surface of neutral spin- 
1/2 fermionic spinous /„ (a =1, |). The physical electron 
operator Ca = fa- Both b and are coupled to a [/(I) 
gauge field a^^. Near the chemical potential tuned Mott 
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transition, the slave boson condensation scale Tf, is of or- 
der the ground state boson phase stiffness pbs- Naively 
once h condenses, by the Anderson-Higgs mechanism the 
gauge fields are gapped out, and the acquire electric 
charge to become Landau quasiparticles. The key point 
of Ref. [H is that, in contrast to this naive expectation, 
the Anderson 'plasnionization' of the gauge field occurs 
only at a scale ~ p^^^ which is much smaller than the 
'Higgs' scale p^s of the b condensation. This is because 
of the unusual dynamics of the gauge field due to Lan- 
dau damping by the fermions. At energy scales above 
the Anderson scale but below the Higgs scale the bosons 
are condensed but the gauge field can still be treated as 
gapless. This leads to strong scattering of the fermions 
and to non-fermi liquid physics. 

We begin by deriving an effective field theory for en- 
ergy scales below the scale T;,. In the absence of spinon 
pairing the full slave boson action (in imaginary time t) 
reads 



becomes 



S 



S[b, a] 



S[f,,,a] + S[b,a] + S[a] (1) 

dr ^ fra [dr 

i(lQ flf) froL (2) 

r 

ts (/™/,.'ce^"'-' + h.c) (3) 

<rr'> 

[dr - iao ~ p.b)br (4) 
tc {brbr'e^"--' + h.c) + V[brbr] (5) 



Here r, r' are the sites of the square lattice. Pb^Pf are 
the boson and fermion chemical potentials respectively. 
They are adjusted so that the condition < brbr >= x, < 
frfr >— I ~ X are both satisfied. ts,tc are the hopping 
amplitudes of the fermion and boson respectively. The V 
term represents a short range boson repulsion. We have 
not explicitly written out the gauge action S[a\ generated 
by integrating out high energy (b.f) modes. For energy 
scales well below Tf, the boson has a well defined ampli- 
tude, and we may write br ~ 5oe'^'' with &o ~ Vx. The 
phase 9r is 'Higgsed' and is described by the Gaussian 
action 

(6) 

where Kb, Pbs ~ x are the compressibility and phase stiff- 
ness of the boson. The index i refers to the two spatial 
directions, and Vi is a lattice derivative. Now we shift 
^ 0.^1 — 9^9 and introduce the fermion field 



-iBr f 

^ J ra 



(7) 



SiFL[q,a\ = 



/ X! ^ ^ ^) 

Jt ^ 

j Yl (gragr'ae'"'-'-' + h.c) (9) 



<rr'> 

2 



E'^b 2 , Pbs -J2 



(10) 



The operator q^a is gauge invariant and has the same 
quantum numbers as the electron operator. Specifically 
Cra = boqra- Thus the truc electron spectral function dif- 
fers from that of the g-field by an overall factor of b^ ~ x. 
The action Si fl i] represents electron-like quasiparti- 
cles coupled to a massive gauge field. It can be under- 
stood as a theory of electrons coupled to strong density- 
density and current-current interactions, as is readily 
seen by integrating out the gauge field. Thus despite 
our derivation using spin-charge separated variables the 
IFL regime is not spin-charge separated. As the coeffi- 
cients Kb, Pbs ~ X the resultant four fermion interactions 
have strength of order \/x and hence are strong. 

We proceed with a standard RPA treatment to deter- 
mine the effective gauge action by integrating out the 
qa- The ao component receives a large x- independent 
constant contribution coming from the bulk compress- 
ibility of the fermions which represents the usual Debye 
screening. A common approximation is to drop the terms 
involving oq, and focus on the transverse component a^. 
It is however more correct to retain it as an effective x- 
dependent on-site Hubbard repulsion Ueff{x). For not 
too small doping this repulsion will not change the qual- 
itative physics but will play a role in determining things 
like the superconducting Tc. The more dangerous trans- 
verse gauge field has an effective action 



S, 



-y 



{\l{p,iujn) + Pbs)\a{p,iujn)\^ (11) 



H, the current-current polarizability of the g-fermions, is 



given (after continuing to real frequencies itOn 
by 



n(p,w) = iuj(T{p,T) ' 
<7{k,T) 



XP 



(12) 
(13) 



An effective action, appropriate below the scale Tb, then 



Here a(0, T) may be interpreted as the conductivity of 
the g-fermions, and 74 is an appropriate (transport) scat- 
tering rate, is the Fermi velocity and the coefficient 
fco ^ Kf (the Fermi wave vector). The diamagnetic sus- 
ceptibility X ~ ^gd^ ~ ^ is of order the inverse fermion 
mass [d is the lattice spacing) . F or T <C the scatter- 
ing rate ^tr found below is sufficiently small that it can be 
ignored compared to the typical momentum p ~ \J 
Then the typical frequency of a gauge fluctuation of mo- 
mentum p scales as Lo . The transverse gauge field 
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a thus gets a gap of order 



coh 



(14) 



where Ep is the Fermi energy. This gap, which is para- 
metrically smaller than ~ pts , sets a scale above which 
the gauge field has still not realized that the boson has 
condensed. 

At intermediate temperatures 3> T 3> Tcoh the q- 
field is strongly scattered by the 'gapless' gauge fluctu- 
ations. It is only on cooling below Tcoh does the gauge 
field become truly gapped, and the scattering subsides to 
produce a Landau fermi liquid. Thus the coherence scale 
is given by Eqn. [TJ] and not by ~ a; as previously 
assumed. 

Despite being metallic, IFL has the same spin physics 
as the insulating spin liquid with a spinon Fermi sur- 
face. Thus it will have a constant spin susceptibility, 
and enhanced "2Kf" spin correlations^ {i.e at wavevec- 
tors connecting tangential portions of the Fermi surface) . 
Further at low temperature (but still above Tcoh), the 
electronic specific heat will have a T^/^ temperature de- 
pendence exactly as in the corresponding insulating spin 
liquid. Interestingly as the Qa field carries both spin and 
charge the enhanced 2Kp correlations will also show up 
in the charge response. 

Consider the single particle Greens function. The q- 
field acquires a non-trivial self energy due to scattering 
off thermally excited low frequency gauge fluctuations. 
In leading order perturbation theory, the scattering rate 
7 at a Fermi surface point K is determined from the 
imaginary part of the q self energy I]"(K, w = 0) and is 
given by 

1 ^ IT [ (VF X p)2 d" (p, ^l)A (K - p, -r!) 
Jp.n 

in{Q) + f{Q)) 

Here vf is the Fermi velocity, D" is the spectral den- 
sity of the gauge field propagator, and A is the spec- 
tral function of the g-field. n(il), f{fl) are the Bose and 
Fermi fimctions respectively. We use A (K — p, — ft) — 
S(—Q — CK-p) where tp is the dispersion of the g-field. 
For K at the Fermi surface, ck-p may be expanded as 
—vfP\\ + cp]_ with pj| , p± the two components of p paral- 
lel and perpendicular to the normal to the Fermi surface. 
The constant c is of order the inverse quasiparticle mass. 
This fermion dispersion implies that in the self energy 
integral the typical jpy | '^-^ «C b_L|- Thus we may drop 
P|| in the gauge propagator and do the integral to get 
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a \sinh{pfl)J Jq 



dpD {p,n) (15) 



The gauge propagator is obtained from Eqns. fTT]andfT2l 
and replacing the polarizability by its zero temperature 
form (justified as the scale for the momentum p is the 
inverse gauge field screening length ^/2mphs ^ ^tr/vp)- 



We may then write the result of the p-integration in the 
following form 



7 = 



where 



git) 



dx- 



n 



sinh{Pfl) 



n 



t^ + {x^ + x) 



(16) 



(17) 



This has the limiting behaviors g{t 0) ~ In j, g{t 
oo) 4-- In the IFL regime we have T > Tcoh- The fl 
integral can be done in this limit and we find a linear T 
single particle scattering rate 



7 - (kBT) 




(18) 



As the scattering off the a fiuctuations is small angle 
the transport scattering rate is potentially different from 
the single particle one. Accordingly Eqn. [15] must be 



modified by an additional factor of 



We then find 



a transport scattering rate jtr ~ T^ (and hence a Ta 
resistivity) at the lowest temperatures. The same rate 
was founc^ for spinon transport in the NFL regime but 
is now interpreted as the transport rate of electrons in 
IFL. 

Next we consider the underdoped region. The mean 
field calculation^ demonstrates that below a tempera- 
ture scale T* a non-zero d^2_y2 pair amplitude develops. 
Within the mean field theory the superconducting transi- 
tion happens at a Tc = min{Tb, Ta)- Beyond mean field 
theory, Tc will be suppressed by various fiuctuation ef- 
fects, as we discuss at the end of the paper. Thus for 
moderate underdoping the pseudogap region (PG of Fig. 
[TJj) must be described as IFL modified by the presence of 
a local d-wave pair amplitude but no global phase coher- 
ence. We now develop a description of the single particle 
properties in the PG regime. 

In the PG regime the action Eqn. [H] must be supple- 
mented by the coupling to a fiuctuating d-wave pair field 
A(r, t). We write 

S[q,A,a] = SiFL + S[A,q] + S[A,a] (19) 
S[A, q] = / X! ^* X! (irTQr'i - g,.ig,.'T)(20) 

+C.C (21) 

Here ry^r' = +1 on horizontal bonds and —1 on vertical 
bonds. We have thus taken A to couple to a c?-wave sin- 
glet pair of the g-fields that has center of mass at the site 
r. We have not explicitly written out the term S[A, a] 
which describes the action for the boson field A coupled 
minimally with charge 2 to the a-field. In a finite temper- 
ature regime in which A has a well-defined amplitude Aq 
but fluctuating phase its correlators will have the form 



{A*ir,t)A{0,0))=AlFir,t) 



(22) 
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such that F(0,0) = 1 while for \t\ > or |r| > 
vp/^, F{r,t) decays exponentially to zero. T is thus the 
Cooper pair phase decay rate. 

Previous work^iSii^ showed how the combination of a 
linear -T scattering rate together with the local d-wave 
pair amplitude lead directly to an electron spectral func- 
tion that has an antinodal pseudogap coexisting with 
gapless T-dependent Fermi arcs near the nodal direction. 
Indeed if Eqn. [TS] survives the development of the pair 
amplitude then the electron self energy has an extra con- 
tribution 



Sp(K,c^) 



A2 



(23) 



which was used to successfully fit the ARPES data^. 

The key microscopic question therefore is to justify 
Eqn. [THl even in the presence of a local pair amplitude. 
We now show that the assumption of such a scattering 
rate is self-consistent within the slave particle gauge the- 
ory. To understand this we first need the form of the 
gauge propagator in this pseudogap region. The polariz- 
ability in Eqn. 1121 will now get contributions from both 
the q and A fields. These may be described in terms of 
a frequency dependent function 



K{uj) L{T)koVF 



(24) 



The first term is the pair field contribution with a 
frequency dependent phase stiffness K{Ld) satisfying 
K{ijj) K, for 1^1 ^ r and K{uj) « pfs (the 'bare' phase 
stiffness of the A field) for P < |w| < Aq. The second 
term in Eqn. [24] is the "quasiparticle" contribution com- 
ing from the gapless Fermi arc region. Its coefficient is 
reduced by a factor L{T) ~ ^ which reflects the frac- 
tional length of the gapl ess arc (in agreement with the 
experiments of Ref. llOl ). For moderately underdoped 
samples L{T) is a sizeable fraction even close to Tc- 

The scattering rate of the q-field at a Fermi point in 
the arc region off these gauge fluctuation is given by Eqn. 
1151 but with modified gauge and electron spectral func- 
tions. Recognizing that the electron spectrum in the arc 
region 7 ^ |Aok| is little modified from that in IFL, 
the p integrals may be done as before. The important 
contribution again comes from |f2| <^T so that 



7 ~ TvpkoL 



{pbs + K{n)) 



■9 



Ep 



(25) 

The function g has most of its support in a region where 
its argument is < o(l). Thus so long as P ^ 

2^ the important frequency range has K{Q) = 0. Then 
the scattering rate continues to be given by Eqn. [181 If T 
falls below this scale the scattering rate will drop as the 
relevant low frequency gauge fiuctuations become stiffer 
due to the contribution from the pair field. 



As T is increased from Tc, we expect that P drops 
rapidly (at least in the critical region within Kosterlitz- 
Thouless theory). Indeed the experiments of Ref. [ill 
show that P ~ T already about 20K above Tc- (How- 
ever P ~ Aq only at much higher temperatures.) Con- 
sequently for moderately underdoped samples the condi- 
tion for incoherent single particle scattering will be sat- 
isfied all the way down to within a narrow temperature 
range of the superconducting transition. Thus our orig- 
inal assumption of a large 7 is self-consistent almost all 
the way upto Tc, and this implies the the emergence of 
both the pseudogap and the gapless T-dependent Fermi 
arcs. 

We emphasize that this mechanism for the incoherent 
scattering in the normal state is distinct from previous 
slave particle theories where the Tc was tied to the bose 
condensation of the holons. Rather here it is the freeing 
up of the low frequency gauge fluctuations by crossing 
Tc that leads to the incoherent scattering. Clearly in 
our theory coherent quasiparticles always exist at T = 
so that a magnetic field sufficient to suppress Tc is not 
sufficient to suppress the quasiparticle coherence scale, 
as assumed in Ref. 

Finally we comment on the effects of the on-site ef- 



fective Hubbard repulsion U, 



effix) 



which should increase 



with decreasing x. First this will lead to short ranged an- 
tiferromagnetic spin correlations, and eventually to anti- 
ferromagnctism for very small x. Second, Ueff will lead 
to a suppression of the superfluid stiffness psf of the pair 
held A from its BCS value. This in turn will bring down 
the true superconducting transition temperature Tc from 
the pairing scale T*. A different mechanism for the sup- 
pression of Tc is due to the thermal excitation of current 
carrying nodal quasiparticle states proposed in Ref. [H. 
Within the present theory the current carried by such a 
quasiparticle is proportional to x for small x. However 
numerical calculations^ with projected wave functions 
show that for moderate x the current rises rapidly. As 
our theory anyway is designed only for such doping levels 
this may be a viable mechanism. We leave for the future 
a more elaborate theory of these effects. 

In summary we explored issues of single particle coher- 
ence and pairing within the slave boson gauge theory of 
the doped Mott insulator. We studied a new non-fermi 
liquid regime, characterized by a linear-T single particle 
scattering rate, which extends to much lower tempera- 
ture than the boson condensation scale. In the under- 
doped side this state acquires a d-wave pair amplitude 
below a T* line, and eventually becomes a c?-wave super- 
conductor. In the regime between Tc and T*, there is a 
pseudogap in the single particle spectrum which coexists 
with 'gapless' Fermi arcs near the nodal direction that 
shrink as the temperature is reduced. Thus this model 
captures many of the most striking phenomenological as- 
pects of cuprate physics. At present we however do not 
have a natural explanation of the linear resistivity that 
is seen at temperatures above T* . Perhaps this difficulty 
can be overcome in the future. 
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